Charge density wave formation near band degeneracies: a new paradigm towards 

broken symmetry ground states 



H.-M. Eiter, 1 M. Lavagnini, 1 R. Hackl*, 1 E.A. Nowadnick, 2 ' 3 ' 4 A.F. Kemper, 2 ' 3 
T.P. Devereaux, 2 ' 3 J.-H. Chu, 2 ' 3 ' 5 J.G. Analytis, 2 ' 3 ' 5 LR. Fisher, 2 ' 3 ' 5 and L. Degiorgi 6 

1 Walther Meissner Institut, Bayerische Akademie der Wissenschaften, 85748 G arching, Germany 
2 Stanford Institute for Materials and Energy Sciences, 
SLAG National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park CA 94025, USA 
3 Geballe Laboratory for Advanced Materials, Stanford University, Stanford CA 94305, USA 
^Department of Physics, Stanford University, Stanford CA 94305, USA 
5 Department of Applied Physics, Stanford University, Stanford CA 94305, USA 
6 Laboratorium fur Festkorperphysik, ETH - Zurich, 8093 Zurich, Switzerland 
(Dated: August 29, 2012) 

Charge and spin density waves, periodic modulations of the electron and magnetization densities, 
respectively, are among the most abundant and non-trivial low-temperature ordered phases in con- 
densed matter. The ordering direction is widely believed to result from the Fermi surface topology 
However, several recent studies indicate that this common view needs to be supplemented. Here, 
we show in which way an enhanced electron-lattice interaction can contribute to or even determine 
the selection of the ordering vector in the model charge density wave (CDW) system ErTe3. Our 
joint experimental and theoretical study allows us to establish a relation between the selection rules 
of the electronic light scattering spectra and the enhanced electron-phonon coupling in the vicinity 
of band degeneracy points. This new proposal towards CDW formation may be of general relevance 
for driving phase transitions into other broken- symmetry ground states, particularly in multi-band 
systems such as the iron based superconductors. 

PACS numbers: 71.45.Lr, 78.30.-j, 63.20.kd, 64.70.K- 



I. INTRODUCTION 

The common view of charge density wave (CDW) for- 
mation was originally posed by Kohn pQ. Following 
Kohn's reasoning, the tendency towards ordering is par- 
ticularly strong in low dimensions, because the Fermi sur- 
face has parallel parts, referred to as nesting. This leads 
to divergences in the Lindhard susceptibility, determining 
in turn the magnitude and the direction of the ordering 
vector Q [2]. However, several recent publications raise 
the question as to whether nesting alone is sufficient to 
explain the observed ordering direction Q [2H5], partic- 
ularly in dimensions higher than one. 

In order to set the stage for our discussion, we 
introduce the two-dimensional rare-earth tri-tellurides 
(RTes), as shown in Fig.[l]A. Among them the prototyp- 
ical ErTe3 undergoes a first CDW transition at Tcdwi 
= 265 K, followed by a second one at Xcdw2 = 155 K 
and allows a robust access to its intrinsic CDW prop- 
erties [6HT5]. The electronic properties of these layered 
CDW compounds can be modeled by considering a sin- 
gle Te plane (Fig.[T]Z? and Appendix). The two dominant 
hopping terms are V pcr and V pn along and perpendicular 
to the overlapping p x or p z orbitals, respectively, lead- 
ing to slightly warped Fermi surface planes (Fig.[l]C). 
For this band structure, Yao and coworkers [3] studied 
the influence of band- filling and electron-phonon coupling 
strength on the charge ordering and established a strong 
coupling limit for the experimentally observed stripe-like 
CDW state. Johannes and Mazin [5 additionally elab- 
orated on the fact that the landscape of the Lindhard 



susceptibility has comparably pronounced peaks at the 
nesting vector Q*, predicted from the band structure, 
and at the CDW ordering vector Qi, determined experi- 
mentally in RTes [8] (Fig.[T]C and Fig.|6]A), quite similar 
to a recent observation in NbSe2 [4]. The inclusion of 
a small but non-zero hopping term between p x and p zi 
Vxz 7^ 0? lifts the degeneracy at the intersection points 
of the two quasi ID Fermi surfaces, as emphasized in 
Fig.[l]D, and improves the agreement with the experi- 
mental Fermi surface [7 , but hardly influences the com- 
petition between Q* and Qi. Therefore, mechanisms be- 
yond purely electronic ones have been conjectured to play 
an important role in selecting the ordering vector and in 
density wave formation. These include orthorhombicity, 
the tendency towards phase separation and nematicity 
via the Coulomb interaction, strongly momentum depen- 
dent electron-phonon interaction due to peculiarities of 
the band structure, breakdown of the Coulomb screen- 
ing and other competing instabilities, such as magnetism 

[3 HUM!]. 



In this study, we analyze data from Raman exper- 
iments and the related selection rules for ErTe3 and 
demonstrate that the lifting of band degeneracies en- 
hances the light-scattering sensitivity and, concomi- 
tantly, the electron-phonon coupling at ordering vectors 
not necessarily coinciding with those predicted by nest- 
ing alone. We show that symmetry may generally play 
a pivotal role in selecting the proper ordering vector in 
multiband systems. 
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FIG. 1: Real and reciprocal space structure of rare earth 
tri-tellurides. (A) Crystal structure with violet and green 
spheres representing tellurium and rare earth atoms, respec- 
tively. The 3D crystallographic cell is indicated in grey. Note 
that the Te planes are perpendicular to the b axis [6 . (B) Or- 
bital character of the Te bp orbit als in the Te plane. The unit 
cell is indicated by a full line. The band structure near the 
Fermi energy Ey is derived from the Te 5p orbitals. The var- 
ious hopping matrix elements are indicated. (C) Fermi sur- 
face as derived from the p x and p z orbitals alone, (V pa - = 2.99, 
Vpir = —1, V xx = 0.09 eV, and V xz = 0, for details see Ap- 
pendix). Note that only V xz contributes to the hybridization. 
There are two energetically possible orientations for CDW or- 
dering, Qi and Q*, where Q* is the wavevector predicted by 
nesting. (D) Theoretical Fermi surface for V xz = 0.12 eV best 
reproducing the experimental findings 0. Also shown are the 
two experimentally observed orthogonal ordering vectors Qi 
and Q2 parallel to the crystallographic c and a axes, respec- 
tively, corresponding to the CDW transitions at Tcdwi and 

TcDW2- 



II. SAMPLES AND EXPERIMENTAL 
TECHNIQUE 

Well characterized single crystals of ErTe3 were grown 
by slow cooling of a binary melt as described elsewhere [6j 
|5]. ErTe3 is a particularly well ordered system. In the a-c 
plane, the resistivity varies only slowly above Tcdwi and 
is very small in the T = limit [8]. The crystals were 
cleaved before being mounted into the cryostat. 

The imaginary part Rx'( s (w, T) of the Raman response 
is measured for various polarization combinations of in- 
coming and scattered photons (e^ and e s ) referred to as 
aa, cc, and ac using Porto notation. Symbolic represen- 
tations by two arrows in the Te plane are shown along 
with the spectra. Usually, more than one symmetry com- 
ponent is projected out at a given polarization e^,e s . 



The pure symmetries correspond to specific eigenvectors 
in the case of phonons and to separate regions in the 
Brillouin zone for electron- hole excitations [20]. In this 
publication we show predominantly spectra with aa and 
cc polarizations which comprise A\ g and B\ g symmetry 
components in a tetragonal lattice and the A g symmetry 
on an orthorhombic lattice. In either case, the response 
has C2 symmetry. 

For the experiments we made use of a solid state laser 
emitting at 532,3 nm (KLASTECH SCHERZO 300) for 
excitation. The absorbed laser power ranged from 1 to 
2 mW to keep the local heating below 5 K in the 50 x 
100 jim 2 sized focus. The spectra were measured with a 
resolution of 2.5 cm -1 at low energy and 20 cm -1 at 
high energy. The Raman response Rx" is then obtained 
by dividing the measured spectra by the thermal Bose 
factor {1 + n((j, T)} = [1 - e -W^]-i. R ig a constant 
which absorbs experimental factors and takes care of the 
units. 



III. RESULTS 

A. Fluctuation regime above the CDW transition 
temperature 

We first display the low frequency Raman spectra 
above Tcdwi i n Fig. [2] The narrow lines superposed 
on the continuum are the Raman-active phonons of the 
high-temperature phase [11] . Instead of the expected flat 
continuum [21], strongly temperature dependent shoul- 
ders, emerging from the normal metallic response, are 
observed for Tcdwi < T < 300 K in the low-energy part 
of the spectra (Figs. [2]A and B). These excitations have 
similar intensity in aa and cc polarization configurations, 
soften and get stronger upon approaching Tcdwi from 
higher temperatures. Above 300 K the spectra are essen- 
tially temperature independent as expected for a metal 
with an almost constant resistivity [8]. 

In Fig.[2]C and D we show the strongly temperature 
dependent parts of the spectra alone which closely follow 
the prediction of Caprara and coworkers for the exchange 
of fluctuations [22]. A fluctuation regime, suppressing 
Tcdwi below the mean field transition temperature and 
out of which the CDW state emerges, is expected to ex- 
ist at temperatures well above Tcdwi •> particularly in our 
case with a large ratio 2Acdwi AbTcdwi ~ 15 being ap- 
proximately four times the canonical mean field value of 
3.53 [2 j. Signatures of such CDW precursor effects were 
indeed observed by x-ray diffraction [8] and by optical 
(IR) spectroscopy [13j[T5]. Here, we further reveal the 
nature and the symmetry properties of the fluctuations; 
our observations are compatible with the A\ g channel, 
which indicates the survival of the C4 rotational sym- 
metry of the pseudo-tetragonal phase [23] . This excludes 
the presence of precursor effects due to nematic order, as, 
e.g., intensively debated in the cuprates p~9j[24], which 
would appear in B\ g symmetry. 
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FIG. 2: Normal state Raman scattering response of ErTe3 at 
low energies. The imaginary part RxI s (oj,T) of the Raman 
response is shown in panel (^4) for aa- and in panel (B) for cc- 
conflguration with the first and the second label representing 
the polarizations and e s of the incoming and, respectively, 
scattered photons as indicated in the insets of panels (C) 
and (D). The spectra display the presence of a fluctuation- 
induced response for a temperature range of about 30 K above 
Tcdwi- In c and d the fits after Ref. [22] to the fluctuation 
contribution Ai?x" = R)d'(u),T) - Rx"(v, 302 K) are shown. 



B. CDW amplitude mode excitations and 
electron-phonon coupling strength 

Immediately below Tcdwi the amplitude mode (AM) 
of the CDW pops up and gains intensity with decreas- 
ing temperature, shown as peak A in Fig. [3] The AM 
appears in both polarizations, aa and cc, with an inten- 
sity ratio of approximately 2:1. Upon further cooling 
the AM moves to higher energies, couples to phonons 
[14], and gains more than one order of magnitude in 
intensity [T2] . At the lowest temperature it saturates 
at uj am = 71 cm -1 (Fig.[3]C, peak A). The relation 
wam = A/A^2k F [2 j at T = between the energy of 
the AM and the un-renormalized CDW phonon energy 
(^2k F = HO ± 20 cm- 1 , Ref. [Tp) leads to A = 0.4 ±0.1, 
indicative of strong electron-phonon coupling. As shown 
in the inset of Fig.[3]C the second AM (peak B) is fully 
cc polarized and saturates at 38 cm -1 . Another mode 
(peak C in Fig.[3]C) with resolution-limited width ap- 
pears at Q§ ~ 18 cm -1 . The full analysis yields tetrago- 
nal Big symmetry, and the energy corresponds to the 
beat frequency of the two amplitude modes £\eat = 
1/2 1 £1 ami — ^AM2| to within the experimental error. Mi- 
croscopically, the coupling between collective states may 
result from eigenvector mixing or because the two con- 
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FIG. 3: Amplitude modes of ErTe 3 for the two CDW tran- 
sitions. The three panels show comparisons of the aa and 
cc polarized spectra at selected temperatures. Note the dif- 
ferent intensity scales. The amplitude mode (peak A) of the 
first transition displays a 2:1 intensity-ratio between the aa 
and cc polarizations, while the second AM (peak B) is fully 
cc polarized. The inset in (C) shows the second amplitude 
mode (peak B) and beat mode (peak C) at 18 cm -1 on an 
expanded intensity scale. 



densates share common parts of the Fermi surface [25]- 
1271. 



C. Temperature dependence and anisotropy of the 
CDW gap 

Figure [4] depicts the electronic Raman response of 
ErTe3 at various temperatures below Tcdwi- The spec- 
tra at 262 K (Fig.[4]A) are isotropic, rise almost linearly 
between 800 and 3500 cm -1 and finally become fiat. 
Upon lowering the temperature there is a transfer of spec- 
tral weight in the aa spectrum from low to high energies 
(Fig.[4]£ and C). At 8K (Fig.[I]c7), there is a relatively 
weak new structure in the cc spectrum in the range 500- 
1200 cm -1 . The insets in Fig.[I]c7 highlight the tern- 
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perature dependences of the aa and cc polarized spectra 
right above and below Xcdwi and Tcdw2, respectively. 
In either case, spectral weight is progressively suppressed 
below the gap edges and piles up above. As opposed 
to the AMs (Fig. [3]), there is a full anisotropy without 
any leakage between the two orthogonal aa and cc di- 
rections, indicating that the crystal is single domain in 
the probed spot. Single domain areas were already ob- 
served in an earlier ARPES experiment on ErTe3 [28] . 
We identify the edges with twice the maximum gap ener- 
gies of the first and the second CDW, 2A 1 = 2800 cm- 1 
and 2A2 = 800 cm -1 , respectively, in agreement with 
ARPES findings [28]. 



D. Raman selection rules and anisotropic 
electron-lattice coupling 

We now elaborate on the remarkable selection rules 
described in the preceding paragraph (Fig. [4| and relate 
them to hybridization effects of the band structure. In 
addition and more importantly, we demonstrate that the 
selection rules and the strong anisotropy of the electron- 
phonon coupling, which influences the CDW ordering, 
are intimately connected and just two sides of the same 
coin. 

The electronic Raman response including the selection 
rules can be derived directly from the band structure and 
the momentum dependences of the CDW gaps using the 
formalism of Vanylos and Virosztek [29]. The intensity 
of the light scattering for different polarization combi- 
nations (Raman vertices) are mainly determined by the 
curvature of the electronic dispersion, as described in de- 
tail in Ref. [20]. If we neglect hybridization (V xz = 0, 
see Fig. [!]£?), the Raman vertices are almost featureless 
with little highlights in particular regions of the Fermi 
surface, as illustrated in Fig.[5]A-C for the aa, cc and ac 
polarizations. Upon including hybridization, the band 
degeneracy is lifted (Fig. T]D), and the two bands exhibit 
strong curvatures (Fig.JlD and Fig. [7] in Appendix. As a 
result, the vertices become highly focused along the di- 
agonals of the Brillouin zone for parallel polarizations, as 
shown in Fig.[5]D, E and G, H for both bands, because 
of nearly singular band curvature [30] . This focusing en- 
hances the light scattering precisely at the Fermi surface 
points connected by the CDW ordering wavevectors. 

The focusing effect on the electronic spectra can be 
demonstrated directly via the weak-coupling Raman re- 
sponse Xr r [ see Appendix Eq. (|8j)]. We assume that at 
T < Xcdw2 the two perpendicular CDWs with ordering 
vectors Qi and Q2 have fully developed gaps Ai and A2, 
respectively. The Qi and Q2 vectors connect the corners 
of the electron pocket around the T point with the cor- 
ners of the hole pockets (Fig. [I] and Fig.|5|. The response 
for the aa, cc, and ac polarization orientations is shown 
in Fig.[5j/-L. No mixing can be observed in the spectra 
with parallel polarizations (panels J and K). In the ac 
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FIG. 4: Temperature dependence of the high energy Raman 
spectra in ErTe3. (A) At 262 K there is no difference in the 
electronic excitations at high energy for both polarizations. 
The insets sketch the incoming and scattered photon polar- 
izations. (B) At 164 K the anisotropy between the two po- 
larizations is already well resolved. (C) In the limit T — >• 
the electronic gaps with edges at 800 and 2800 cm -1 (dashed 
lines) for both CDWs are fully developed. The upper left and 
lower right insets show the temperature dependences of both 
CDW gaps, using the same colour-code. Note that the sig- 
natures of the CDW transition at Tcdwi with Qi || c are 
observed with incoming and scattered light polarizations par- 
allel to a. The opposite is true for the transition at Tcdw2 as 
shown in detail in Fig. [5] 



configuration both gaps are in principle visible (panel L) , 
but the expected intensity is three orders of magnitude 
smaller than in aa and cc and cannot be observed in the 
experiment. 

As we do not include any other scattering mecha- 
nisms, phase space restricts the non-resonant creation of 
electron-hole excitations to points where the CDW mixes 
particles with wavevectors k and k + Q across the Fermi 
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FIG. 5: Theoretical prediction for the Raman vertices and 
spectra of ErTe3. (A-I) All vertices are derived from the band 
structure according to Eq. Q in the Appendix. (A-C) The 
first row shows the vertices without 2D coupling (V xz = 0) 
and the corresponding Fermi surfaces. (D-F) Raman vertices 
for the bands corresponding to the central part of the Fermi 
surface (green) and (G-I) to the outer part (light blue). The 
focussing effect due to the lifted degeneracy enhances the ver- 
tices by more than two orders of magnitude as indicated by 
the colour code. The ordering vectors are displayed in the 
third row. (J-L) All spectra are calculated at T < Xcdw2 
and include both CDWs (Eq. |8} in the Appendix). (J, K) 
For parallel polarizations one observes only the CDW with 
ordering vector perpendicular to the light polarizations. The 
response of the respective orthogonal CDW is too weak to be 
visible. (L) For ac polarization both gaps can be resolved but 
the overall intensity is more than three orders of magnitude 
lower than that in the two other configurations. This is well 
below the detection limit, and, in fact, no signatures of the 
gaps can be observed experimentally in ac polarization. 



surface. Consequently, light scattering is enhanced where 
energy is gained due to the CDW gap opening at the 
Fermi surface (Fig.[5|/-L). Raman scattering efficiently 
projects out the relevant parts of the Brillouin zone in 
such a multi-band system, so that the signal is small at 
low energies below the gap edge but is significantly en- 
hanced at twice the CDW gap (Fig.[5|/-L), in agreement 
with our experiment (Fig.[4](7). 

Obviously, the lifting of band degeneracies dramati- 
cally affects the Raman selection rules by locally enhanc- 
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FIG. 6: Comparison of the susceptibilities. (A) 3D plot of 
the real part Xl °f the Lindhard susceptibility. There is little 
structure around the rim, but Q* (blue arrow) is the prefer- 
able ordering vector. (B) 2D superposition of the real parts 
of the Lindhard susceptibility (xl, yellow-red) and the pro- 
jected susceptibility for interband scattering transitions (x'p, 
green). It is the focusing effect of the stress tensor which 
selects the experimentally observed ordering wave vector Qi. 



ing the Raman vertex (Fig. [5]) which is proportional to 
the inverse effective mass tensor [20] . Now we make use of 
an analogy between electron-phonon and electron-photon 
scattering, where the electrons scatter from a phonon 
rather than a photon and the electron-phonon coupling 
vertex replaces the Raman vertex [31 . In particular, for 
the case of stress phonons, the electron-phonon coupling 
vertex is given by the electronic stress tensor, which is 
proportional to the inverse effective mass tensor [32H35] . 
As a result, for the system studied here, the electron- 
phonon coupling vertex, like the Raman vertex, is highly 
anisotropic in momentum space. While the Lindhard sus- 
ceptibility xl is relevant for a momentum-independent 
electron-phonon coupling, in this case we must include 
the momentum dependent electron-phonon coupling ver- 
tex into the electronic susceptibility; we call this suscep- 
tibility xp ( see Appendix) [36 j. Fig. [6] illustrates the im- 
portance of including the electron-phonon coupling ver- 
tex, where we compare the real parts x'l anc ^ Xp °f both 
susceptibilities. Whereas the Lindhard susceptibility x'l 
has maxima of comparable height for several different or- 
dering vectors (Fig.[6]A) and therefore does not lead to 
an unambiguous selection of one of them, x'p contributes 
to the instability at the proper location in q space and 
finally selects the experimentally observed ordering vec- 
tor Qi (Fig.[6]£?). Here, q is the difference of the mo- 
menta k and k 7 of a scattered electron. Furthermore, as 
noted by Yao et al. [3 , any enhancement of the aver- 
aged electron-phonon coupling strength A (definition see 
Appendix) will drive the system further towards the ob- 
served order. These two effects conspire to minimize the 
dependence on model details. 



IV. DISCUSSION 

The huge modulation of the Raman vertex indicates 
strong fluctuations of the charge density in the vicinity 
of the degeneracy points. These fluctuations manifest 
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themselves in the Raman response above Tcdwi (Fi, 
where the lattice still has the full C4 symmetry 
Yao and coworkers showed in the framework of a Lan- 
dau expansion of the free energy [3] that both charge 
fluctuations and electron-phonon coupling cooperate to 
drive the system towards the experimental ordering vec- 
tor. The substantial charge fluctuations couple to and 
soften the phonon relevant for the CDW ordering. Tak- 
ing into account the large electron-phonon coupling near 
the band degeneracy points (Fig. [5]), the charge fluctua- 
tions actually renormalize the phonon frequency at the 
momenta q = Qi(Q2) rather than Q*. Below Tcdwi ? 
the weak orthorhombicity along with the relatively large 
electron-phonon coupling [3] then tips the balance, and 
Qi aligns along c rather than a. Finally, upon further 
lowering the temperature below Tcdw2 Q2 aligns along 
a since the Fermi surface along the c direction is already 
fully gapped by the first transition [HI [28] . 

Therefore, we identify two cooperating effects deter- 
mining the overall selection of the ordering vector: (i) 
The system gains energy by gapping the band degener- 
acy points on the Fermi surface where the Raman se- 
lection rules indicate substantial fluctuations with four- 
fold symmetry above Tcdwi- (h) Since the electron- 
phonon coupling vertex is proportional to the Raman 
vertex for stress phonons 32. 34 . both quantities are en- 
hanced near band degeneracies. For small hybridization 
and an electron-phonon coupling strength of A > 0.5 the 
focussing effect may even be the most relevant contribu- 
tion to the phonon renormalization, and thus the CDW 
formation, while it is only a correction for larger V xz . 
Hence, while electron-phonon coupling is known to be 
important in CDW systems [3 , we identify here on a mi- 
croscopic basis the focussing effect to be a more generic 
paradigm for multiband materials. In this context, it 
seems particularly interesting to explore this novel sce- 
nario in materials such as the iron-based superconductors 
[30] . in which density- wave order and superconductivity 
coexist on a microscopic scale. 
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APPENDIX 

1. Raman scattering and band structure 

For the case of non-resonant light scattering, the Ra- 
man cross section in the limit of small momentum trans- 
fer is given in terms of correlation functions of an effective 
charge density [20] . 



^2 7^( k ^^ S )4,k,<7 C ^k,a, 



(1) 



where cj, kcr (c u ^ a ) creates (removes) an electron with 
momentum k and spin a from band v. In principle, the 
scattering amplitude 7^ depends on both the incoming 
(scattered) light polarizations e^ s ^ and frequencies oj^ s y 
However, in the limit where u^( s ) are much smaller than 
any relevant interband transition frequency, the scatter- 
ing amplitude simplifies to the well known effective mass 
tensor, 



0) = h £ e « 



a,fi=x,y,z 



G> 2 e^>(k) 
dk a dkp 



"a- 



(2) 



Thus the curvature of the bands and the light polariza- 
tion orientations determine which carriers are involved 
in light scattering in different bands and regions of the 
Brillouin zone. 

The tri-telluride band structure for the two bands cut- 
ting the Fermi level is well described by a tight binding 
model that only includes the Te p x and p z orbitals [7]. 
All expressions in this paper are given for coordinates in 
the unit cell defined by the square lattice of Te atoms 
(Fig.[l]B). 

The real space Hamiltonian can be expressed as 



m,n(px) + ' C m+l,n(pz) C "m,n(pz) 



: m,n(px) ~7~ C m,n+l(pz) C ™Mpz) 



H — V p(J ^2 [ C ln,n+l(px) C ' 

— Vxx ^2 [ C L+l,n+l(a) C ^^(a)+ C ]n-l,n+l(a) C ' 

m,n,a 

— Vxz^2 \^rn+l,n+l(px) C ™,n(pz) + C ra+l,n+l(pz) C 

m,n 

+ V XZ ^2 \^m-l,n+l(px) C ™MP z )^^rn-l,n+l(px) C ™, 



ra,n(ai) 



m,n(px) 



\1 y ^ °m,n,a°ra,n,a 
m,n,a 



+h.c. 



n(pz) 

(3) 



where , ^ creates an electron at site (m, n) in orbital 
a = Px,Pz- It includes both nearest and next-nearest 
neighbor hopping matrix elements of one tellurium plane, 
as denoted in Fig. [T]R Fourier transforming Eq. ^ leads 

4" r\ r\ 4--ittv~v Konrlo tttV> nrn 7/ I and V 1 o V\ r\ 1 



to the two bands , where v 



label 
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the electron and the hole band, respectively, 



hl 1 + h™±J{h£-h?y + A{h}?y 



(4) 



With the Te-Te distance set to unity, the energies h 1 ^ 
read 



hi 1 = 
hi 2 = 



h 



12 



-2V p7r cos k z — 2V pa cos k x — W xx cos k z cos k x —/ji, 
-2V p7T cos k x — 2V pa cos k z —AV xx cos k z cos k x —fi, 
-4V XZ s'mk z s'mk x . (5) 



Using the experimental distances between rare earth 
and the Te ions of 4.34 and 4.34/>/2A, respectively, 
the hopping matrix elements are V pa — 2.99, V^- = — 1 
eV. The diagonal hopping integrals are chosen to be 
V xx — 0.09, V xz = 0.12 eV. This parameter choice is moti- 
vated in order to most closely match the ARPES-derived 
Fermi surfaces from Ref. [25] . The chemical potential is 
chosen to yield a filling of 1.6 electrons/band. Any small 
interaction (here we focus on inter-orbital hybridization 
Vxz) leads to sharp corners in the electron and hole band, 
displayed separately in Fig. [7]A and 5, respectively. This 
hybridization thus lifts the band degeneracy along the 
r — M Brillouin zone cut with respect to the V xz = 
case, as shown in Fig.[7]C. 

Finally, we shall note that a more accurately derived 
tight-binding representation would require the inclusion 
of further hopping matrix elements that could affect the 
overall curvature of the bands. Furthermore, another 
larger unit cell taking into account the Te double layers 
is usually considered for addressing the ARPES results 
[7J. While these latter issues would allow a fine tuning 
of the calculations, we trust that this would not alter 
the main conclusion of our work, effectively based on a 
simplified tight-binding approach for a small unit cell. 



2. Weak-coupling Raman response 

The Hamiltonian for a system with two CDWs can be 
expressed in Nambu notation as Hqbw — -^(Qi? Ai) + 
iJ(Q 2 , A 2 ), where 



tf(Q,A) = ^^A k (A)Vw 



(6) 



kcT 



This Hamiltonian is defined in terms of the spinor ip^a- = 

(c+,k-Q,cn c +,k,cr, C +5 k + Q 5Cr , C_ 5 k_Q 5Cr , C_ ? k,cr, C_ 5 k + Q 5Cr ) 

where c±^ a destroys an electron in band with 
momentum k and spin a. The matrix Ak is given by 
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FIG. 7: Band structure of ErTe3 including diagonal hop- 
ping. Panels (A) and (B) show the electron and hole band, 
respectively, using V xx = 0.09 eV and a hybridization of 
Vxz = 0.12 eV. Both the 2D band structure and the constant 
energy contours are plotted. The sharp corners result from 
the hybridization V xz and lead to the huge enhancement of 
the Raman vertices shown in Fig.[5]D-/. Panel (C) shows the 
dispersion along the line r— M— X— Tas indicated in the inset 
before (dashed lines) and after (full lines) the hybridization 
Vxz has been turned on. 



where the diagonal elements correspond to the quasipar- 
ticle energies at k — Q, k, and k + Q, and A is the 
CDW gap. For Q and A we use Qi = (2.16, 0, 2.16)/a, 
Q 2 = (2.16,0, -2. 16) /a, and Ai = 1400 cm- 1 and A 2 = 
400 cm -1 , respectively. The expansion is truncated af- 
ter the first harmonic although higher order terms nQ 
should be included for the incommensurability of Qi 2 

0. 

Computing the matrix Green's function G(k, iuj n ) — 
(iuu n — Ak) -1 and taking into account the Raman vertices 
leads to the Raman response [20], [29] 

Xr,r(q = 0,iO n ) = (8) 
2 Tr[f k G(k,^ n )f k G(k,^ n + ^ n )]. 



Here, itt n and iuj n are the external and internal Mat- 
subara frequencies, respectively, j3 = (/c^T) -1 , and Tr 
denotes the trace. 



The Raman vertex f k is 



/7k + - 



r k 



Q 







7k^ 











7kV Q 











7k-Q 











7k 



(9) 



7k +Q / 
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with 7 ± labeling the effective mass vertices [Eq. de- 
rived from e k [Eq. Q]. Eq. ^ can be evaluated by di- 
agonalizing the Green's functions before calculating the 
trace, which simplifies the analytic structure of the ex- 
pression and allows us to express the response as a sum 
of Lindhard-like terms. 



3. Effect of near band degeneracies on the phonon 
spectrum 

We recall that in a conventional electron-phonon cou- 
pling scenario, a divergence in the electronic susceptibil- 
ity is conveyed via the electron-phonon coupling into the 
phonon dispersion, a phenomenon known as the Kohn 
anomaly pQ. With increasing interaction A, the phonon 
frequency cj q of a branch is first renormalized in a wide 
range of momenta, then softens and finally exhibits a dis- 
continuous derivative at wave vectors Q that nest the 
Fermi surface. A phonon that softens to zero frequency 
at a certain wave vector thus signifies the formation of 
a static lattice distortion, or CDW with that ordering 
vector. 

The renormalized phonon frequency u) q is given by [2] 



where 



A(q, w q ) 



;[1-A(q,w q )], 



25 x'(q,w q ) 



(10) 



(11) 



x'(q, cj q ) is the electronic susceptibility, g in Eq. (11) 
sets the overall strength of the electron-phonon interac- 
tion. The phonon renormalization is expected to be quite 
relevant in the tri-tellurides, since Raman experiments 
have established a strong connection between the CDW 
gap and lattice distortions [TT] . 

We now examine the effect of charge fluctuations on 
the phonon renormalization. We start with the original 
proposal by Kohn, pQ where the Lindhard susceptibility 
Xl directly couples to the phonon propagator with a bare 
charge vertex, 



Xi (q,tt) = 2^ 



/(*£+,) 



f(4) 



Q + iS + 6 



k+q 



J' 



(12) 



Here, a and /3 label the bands. The effect of the normal 
charge susceptibility xl has been studied in some detail 
by several authors [3j[5]. However, within this picture, 
the selection of the ordering vector is not pronounced; the 
maxima in the susceptibility are almost equally high for 
a large manifold of ordering vectors (Fig.[6]A), support- 
ing the conclusion that another mechanism is responsible 
for finally selecting the experimentally observed ordering 
vector Qi. 

We thus consider a generic coupling to stress phonons 
where the electron-phonon coupling vertex is determined 
by the inverse effective mass tensor, which is related to 




FIG. 8: Renormalized phonon dispersion relation. (A) Renor- 
malized phonon frequency (cj q /cj q ) 2 with electron-phonon 
coupling A q =o = 0.5 and diagonal hopping terms set for clar- 
ity to 10 times smaller values than those used in Figs.[T]D and 
[5] The arrow represents the ordering wavevector Qi as de- 
termined experimentally. (B-D), Cuts along the line q x = q z 
(dashed line in panel A) for various values of the electron- 
phonon coupling starting at A = 0.01 (top curve) and then 
ranging between A = 0.05 and A = 0.65 (bottom curve) with 
incremental steps of 0.1. Each panel is for a different value 
of the hybridization parameters: (B) V xx /10, V xz /10; (C) 
V xx /4, V xz /4; and (D) V xx , V xz . 



the Raman vertex via g^ = g^ [32j[35]. As a result, in re- 
gions of strong band curvature the electron-phonon cou- 
pling is large. We consider only interband contributions 
which are generally accentuated by nesting as shown in 
Fig.[l]C. Using the mass density operator [Eq. 0] we de- 
fine the projected susceptibility xp f° r inter-band scat- 
tering transitions as 



k 



/( e k+q/2) /( 6 k-a/2) 



z^7 k 7k 



k-q/2; 



Q + iS + e 



(13) 



k+q/2 k-q/2 



where r± 



k are the two bands defined in Eq. 



4b and 



7k = 7aa + 7cc are ^ ne fr^Y symmetric_effective mass 



vertices derived from these bands [Eq. ([2|]. Note that 
this quantity is only C4 symmetric in the larger momen- 
tum range q^) = (0, 47r) and vanishes for V xz = 0. We 
use this susceptibility in our calculation of the renormal- 
ized phonon frequency [Eqs. ( 10 ) and ( 11 )] . Here we have 
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taken the bare phonon frequency cj q to be independent 
of momentum, hence cj q = cj2q F , as any mild momentum 
dependence will not affect the general behaviour near Qi. 

The phonon frequency renormalized by %p(q, fi) is 
shown in Fig.[8]A for small hybridization (i.e., with values 
of V xx and V xz 10 times smaller than those employed in 
Figs. [i]D and [5| in order to emphasize the effect. Pan- 
els B-D further display the phonon softening along the 
(0,0) — (7r,7r) Brillouin zone cut (dashed line in panel 
i) as a function of V xz and A. It is evident that the 
phonon frequency exhibits a significant softening at the 
wavevector Qi connecting the corners of the Fermi sur- 
face, where the band degeneracy is lifted by hybridiza- 



tion. This predicted phonon renormalization remains to 
be experimentally confirmed by neutron or inelastic x- 
ray scattering (RIXS), since Raman scattering is limited 
to small momentum transfer q <C Qi. The focusing ef- 
fect of the large curvature of the energy bands given by 
the mass tensor is very sensitive to the hybridization pa- 
rameter V xz (Fig.^p-D). Nonetheless, our study has a 
model character pointing out the relevance of the pro- 
jected susceptibility for interband scattering xp- The 
exact strength of the enhancement due to the xp a t the 
experimentally observed ordering vector Qi depends crit- 
ically on the dispersion over the whole Brillouin zone [5] 
which is left to a future study. 
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